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Abstract 

Lattice mismatch in epitaxial layered heterostructures with small characteristic lengths 
induces large, spatially nonuniform strains. The components of the strain tensor have been 
shown experimentally to affect the electronic properties of semiconductor structures. Here a 
technique is presented for calculating the influence of strain on electronic properties. First the 
linear elastic strain in a quantum dot or wire is determined by a finite element calculation. 
A strain-induced potential field that shifts and couples the valence subbands in the structure 
is then determined from deformation potential theory. The time-independent Schrodinger 
equation, including the nonuniform strain induced potential and a potential due to the het- 
erostructure layers, is then solved, also by means of the finite element method. The solution 
consists of the wavefunctions and energies of states confined to the active region of the struc- 
ture; these are the features which govern the electronic and transport properties of devices. As 
examples, two Sia;Gei_2; submicron resonant tunneling devices, a quantum wire with two di- 
mensional confinement and a quantum dot with three dimensional confinement, are analyzed. 
Experimentally measured resonant tunneling current peaks corresponding to the valence sub- 
bands in the material are modeled by generating densities of confined states in the structures. 
Size and composition dependent strain effects are examined for both devices. In both the 
quantum dot and the quantum wire, the strain effects on the wavefunctions and energies of 
confined states are evident in the calculated densities of confined states in the structures, which 
are found to be consistent with experimentally measured tunneling current/ voltage curves for 
resonant tunneling diodes. 
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I. INTRODUCTION 



Epitaxially grown semiconductor heterostructures often consist of materials with lattice param- 
eters that are mismatched by as much as several percent. For thin films of large lateral extent, 
these strains are spatially uniform and the effects are well understood. However, structures of 
relatively small lateral extent, having distinctive geometric features and bounded by free surfaces, 
are characterized by strains that are highly nonuniform. The effects of nonuniform strain on the 
electronic properties of semiconductor heterostructures have been observed experimentally, but 
the coupled physical phenomenon has not been extensively modeled. The analysis of strain 
effects in a quantum mechanical model of semiconductor devices has only recently been attempted 
by Pryor et al. and Williamson et al. ^, who calculate strain induced potentials and wavefunc- 
tions in quantum dots, and Zunger^, who reviews the topic of electronic structure in pyramidal 
semiconductor quantum dots based on atomistic methods. 

Strained semiconductor devices that are based on quantum effects, particularly charge con- 
finement in one or more spatial dimensions, underlie a potentially significant technology. Much 
can be learned about quantum effects by studying the class of devices based on resonant tun- 
neling of carriers from an emitter region, into a quantum well, and then into a collector region. 
The material combinations in these devices and their geometrical features, including the layered 
structure and the free surfaces, lead to complicated mechanical strain fields. Because the sequen- 
tial tunneling of carriers is a simple phenomenon governed by the spectrum of available states in 
the quantum well, and because the devices are extremely small and the operating temperatures 
are low, it is likely that the effects of strain on the electronic and transport properties can be 
represented quantitatively through modeling. 

Calculations of elastic strain fields in semiconductor structures are well suited for the finite 
element method, which is a common tool in continuum mechanics.^'^° In the technique presented 
here, the finite element calculation of the strain in a device is made using a general purpose 
structural mechanics finite element package. However, the use of the finite element method 
(FEM) in quantum mechanics, which is reviewed by Linderberg is much less common. Mod- 
els of semiconductor devices by means of FEM have been proposed by a number of authors for 
one and two dimensional problems. Several FEM models are available for one dimensional reso- 
nant tunneling structures which include the effects of arbitrary potential profiles due to layered 
composition. Chen^*^ models a one dimensional resonant tunneling diode using FEM and 
calculates a current-voltage curve based on a quantum hydrodynamic model. Electron wave- 
functions and band structures for two dimensional quantum wires or quantum dots are analyzed 
using FEM by a number of investigators^ ''^^^ . Tsuchida and Tsukada^^ calculate the electronic 
structure of a perfect Si lattice using a three dimensional FEM formulation. However, none of 
the studies consider a three dimensional device, and none have considered the effects of strain on 
experimentally measured electronic properties of devices. 

The finite element method is well suited for finding approximate solutions of boTindary value 
problems for partial differential equations in finite domains, especially if the (unknown) exact 
solution is a minimizer of a total energy functional. Both the stress boundary value problem 
and the quantum mechanical boundary value problem in the present study are of this type. The 
central idea of the method is that an unknown continuous field in the domain is represented 
approximately in terms of its values at discrete points (nodes) within the domain; the goal is 
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to determine optimal values for these nodal quantities. The domain is covered with areas (in 
two dimensions) or volumes (in three dimensions) whose boundaries are defined geometrically by 
the nodes; these areas or volumes are the elements. Fields are defined within each element in 
terms of the values of the nodal quantities on its boundary by means of a suitable interpolation 
scheme. With the complete field defined in terms of the nodal quantities, the total energy can be 
expressed in terms of the global vector of nodal quantities. Imposition of the condition that the 
actual values of the nodal quantities must render the total energy a minimum leads to a system 
of algebraic equations for these values. The method is ideally suited for numerical analysis by 
computer. In general, it is convergent as nodal spacing diminishes for elliptic partial differential 
equations which have unique solutions. It should be noted that the method is of far broader 
applicability than is implied by these introductory comments. 

In this work, a finite element model is used to analyze both the continuum mechanics and 
the quantum mechanics of a strained semiconductor device. The strain field is shown to affect 
the performance of the device. The model is used to analyze devices studied experimentally 
by Akyiiz et al. ^ and Lukey et al. ^, who show the effects of nonuniform strain on the resonant 
tunneling current vs. voltage characteristics of Si^Gei-a; quantum wires and quantum dots, shown 
schematically in Figure 1. The main features of the model are described in Section II; these include 
the strain calculation, the treatment of the strain effects by deformation potential theory, and 
the quantum mechanical model. The finite element formulation for the two or three dimensional 
Schrodinger equation incorporating valence subband coupling and the nonuniform potential is 
outlined in Section III. The results of the calculations for the resonant tunneling structures are 
discussed in Section IV, and comparisons are made between the calculations and the available 
experimental results. 

II. CONTINUUM AND QUANTUM MECHANICAL MODELS 

The analysis of a strained semiconductor hetcrostructurc is divided into three caknilations. First, 
a linear elastic finite clement calculation is made to determine the strain field, which is a function 
of the composition and the geometry of the structure. Second, the strain induced potential field is 
calculated using deformation potential theory. Third, the time-independent Schrodinger equation 
including the strain induced potential is solved numerically by means of the finite element method 
to obtain the spectrum of energies and wavefunctions of available states. 

A. Strain Field Calculation 

The strain field due to the constraint of epitaxy associated with the mismatched lattice parameters 
of the heterostructure layers is determined within the framework of linear elasticity theory. The 
structure is discretized spatially with a mesh which is more refined near free surfaces and in regions 
where the mismatch between adjacent layers is larger. The finite element mesh used for a strain 
calculation in the resonant tunneling diode quantum wire is shown in Figure 2. The mismatch 
condition is imposed in the finite element calculation by prescribing in each layer a uniform stress- 
free dilatation that is proportional to the bulk lattice parameter of the material in that layer. 
Continuity of displacements is required across the layer interfaces; it is this constraint which 
gives rise to stress. The outer surfaces are considered to be traction free, implying that certain 
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components of the stress tensor vanish, and the material is allowed to relax until it reaches an 
equilibrium configuration. Thus, using a standard structural mechanics finite element program 
the complete state of stress, strain and displacement is determined throughout the device. Strain 
components are shown in Figures 3 and 4 for the resonant tunneling diode quantum wire. Strain 
is a tensor quantity, and its components show significant variation with position throughout the 
device. In smaller devices, the strain is more nonuniform due to the pervasive effect of relaxation 
at the free surfaces. 



B. Strain Induced Potential 

The components of strain induce a potential field that affects the wavefunctions and energies 

of the charge carriers in an otherwise perfect crystal. From first order perturbation theory, the 
strain induced potential that affects wavefunctions in subbands a and /? is formed from the tensor 
product (Singh^^) 

VrP{f)=D,j-f'eij{f) (1) 

where r is an arbitrary position vector, eij{r) is the strain tensor field, and D°j^ is the deformation 
potential tensor for subbands a and /?, which consists of components derived experimentally. The 
indices ij range over the coordinate directions. For the SiajGei-^: material combination, the 
valence band electronic properties are dominated by the heavy hole and light hole subbands, so 
the aP basis includes the heavy hole (HH) subbands denoted by | |, ±| >, and the light hole (LH) 
subbands denoted by | |, ±^ >• The split-off subbands (SO) are ignored because the separation 
energy is considered to be large enough so that coupling effects can be neglected. The potential 
equation can be written as 
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where the terms Dfj^ are the deformation potential tensors which range over the spatial dimensions 
i and j. The deformation potential tensors contain material constants which vary spatially as 
the composition varies in the device; values of these constants for a wide range of materials are 
known from experiments. The repeated ij indices indicate the scalar product (contraction over 
the product of second rank tensors) between the deformation potential tensor and the strain 
tensor. Details of the D'^j^ terms are given in Appendix A, along with the material constants used 
in the calculations. Thus, for a calculated strain tensor function of position ey (r), it is possible 
to calculate the deformation potential function of position V^^ (r) to be included in the quantum 
mechanical analysis. 



C. Quantum Mechanical Model 

In resonant tunneling structures, tunneling currents are determined by available quantized states 
for individual charge carriers. The energies and wavefunctions of a single carrier in the semicon- 
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ductor structure are solutions of the time independent Schrodinger equation 



Hfp{r)^^{r) + y"/^(f)^^(f) = £;^"(f) (3) 

where ^'"(r) is the wavefunction in subband a, E is the energy, H'j^l^{r) is the k • p Hamiltonian 
operator, and V^"^(r) is a potential function of position. 

The k • p perturbation method is used to model the medium. Using this technique, there 
is a tensor function for the effective mass associated with each subband, and there are k • p 
terms coupling the effective masses in different subbands. Written in the same form as dH) the 
Hamiltonian is 
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where the tensors on the diagonal of the matrix are the effective mass tensors for each 

subband, and the off diagonal L^^ tensors introduce a k • p coupling of subbands. Like the terms 

in the deformation potential tensors, the components of L^^ contain material constants and thus 

vary spatially throughout the device. The exact forms of the L^^ tensors, which come from the 
Luttinger-Kohn Hamiltonian, are given in Appendix A. 

The deformation potential and effective mass material constants in @ and are functions 
of the local composition of the device. Values for the constants in each layer are given by linear 
interpolation of the material constants associated with each of the pure elements composing that 
layer. 

The nonuniform potential y^^(r) consists of contributions from two sources and is given by 

^ yc.fi ^ yap 

where V^"^(r) is due to the valence band alignment of material at a given position in the device, 
and V^^{r) is the strain induced potential given in (||). The strain induced potential y^"^(r), like 
the components of the strain tensor, is in general nonuniform in both the lateral directions and 
the vertical direction in the structure. The composition based band offset potential V^^{r) is 
nonuniform only in the growth, or vertical direction in the structure. The total potential y°^(r) 
is shown in Figure 5 for a representative quantum dot calculation. 



III. FINITE ELEMENT TECHNIQUE FOR THE SCHRODINGER 
EQUATION 

A. Finite Element Formulation 

The physical domain of the device is discretized into a mesh of nodes and elements. Elements 
used for the quantum wires and quantum dots are described in Appendix B. The mesh extends 
to the free or insulated surfaces, which impose an infinite potential on the wavefunction. The 
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wavefunctions and energies of the states localized in the active region are insensitive to remote 
boundary conditions, i.e. conditions at boundaries located a large distance away relative to the 
active region size. The mesh is more refined in the active region of the device, where large 
wavefunction gradients are expected. 

The form of the Schrodinger equation to be solved on the finite element mesh is obtained by 
minimizing the total variation of the weak form of the equation with respect to the wavefunction. 
The minimum in variation with respect to the wavefunction corresponds physically to a mini- 
mum energy. Details of the complete variational formulation of the finite element technique are 
included in Appendix B. The functional corresponding to the weak form of the time independent 
Schrodinger equation with a nonuniform potential is given by 

n(^'") = / V^'°L°^V*'^(ii2 + / ^"y^^^^di? -E f ^''^'^dR (6) 
Jr Jr Jr 

where L"^, and V"'^ are functions of position in the structure. The term L"^ is taken to 
be constant within each element of the mesh. The fields V^'", and V"'^ are represented by 
their nodal values. Values throughout each element are determined by interpolation according to 
the particular shape functions that are adopted. Thus, for shape functions N(r) used here, these 
fields are written as 

all nodes all nodes all nodes 

A=l B=l C=l 

where the coefficients in the summations are the values of the fields at the individual nodes. The 
functional n(\I'") is rewritten in a discrete form using (^) and then minimized with respect to 
the nodal values of the wavefunction, Integrals over the region R of the entire structure are 
replaced with integrals over individual element volumes (O*^) and a summation over all elements, 
so that the finite element form of the Schrodinger equation becomes 



elems 

E 



elems 



E E / "^NaL^^vnb + NaNb Y vfNc ]dn ^'1 = eY E E / NANBdn 

(8) 

which is a form ideally suited for computation. 

This form of the Schrodinger equation can be much more simply expressed as 

Kij^j = EMij^j (9) 

where the indices i and j range over all nodal wavefunction degrees of freedom, and the repeated 
indices indicate a summation. Details of the assembly of the matrices K and M are given in 
Appendix B. 
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B. Finite Element Solution 

1. Energies and wavefunctions 

The finite element expression of the Schrodinger equation (|9|) is in the form of a generalized 
eigenvalue equation. The problem has na solutions, where n is the total number of nodes in the 
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mesh and a is the number of subbands in the quantum mechanical basis. The solutions consist 
of energies E and wavefunctions of which the lowest energy states are the most accurate. 

Of the na eigenstates, some states can be found for which the wavefunction ^'(r) or the 
probability density |^'(r)p is confined to the active region of the device. Examples of eigenstates 
confined to the quantum well layers of the quantum dot are shown in Figure 6. The lowest energy 
states of these confined eigenstates are the most relevant to transport. For the example of the 
resonant tunneling diode, sequential tunneling through the double barriers is possible only when 
the tunneling carriers have energies equal to the energies of the confined states in the quantum 
well. Thus, over a range of applied biases, the excited carriers can access the confined states and 
induce a tunneling current only at certain resonances corresponding to the spectrum of eigenvalues 
given by the finite element solution. 

The energy and wavefunction solutions reflect the effects of strain, composition, and effective 
mass on the carriers. The valence band offset in adjacent layers imposes a large relative potential 
on the charge carriers which results in confinement to the quantum well region of the device. The 
strain induced potential is considerably smaller than the valence band potential, but it also shifts 
the wavefunctions spatially and energetically. Confined states corresponding to valence subbands 
with higher effective masses occur at lower energies. 

2. Density of states 

The density, with respect to energy, of states confined to the active region of the device can 
be obtained directly from the spectrum of eigenstates given in the finite element solution. This 
density of confined states is a real space measure of the electronic properties of the device. Effects 
due to strain, composition, size, and device characteristics can be seen in the density of confined 
states. For the example of the resonant tunneling diode, the density of confined states can be 
used to examine strain effects on the resonant tunneling spectrum. A larger density of confined 
states with a given energy increases the probability of resonant tunneling by carriers with the 
same energy. A high probability of tunneling at a given energy is measured experimentally as a 
tunneling current peak. Thus, the density of confined states can be used to make contact with 
experimental data. 

A Gaussian broadening technique can be used to calculate a continuous density of confined 
states function p{E). The energy Ei of each confined state is broadened by a narrow, normalized 
Gaussian distribution, and a summation over all n states gives p{E) as 



where a is a free parameter that controls the width of the Gaussian distributions and thus the 
smoothness of the density of states p(E). The parameter a is chosen to bring out the general 
features of p{E); the broadening of each state is larger than the typical separation of individual 
eigenstates, but narrow enough to bring out features of the density of states that are due to small 
groups of related eigenstates. Typical values of a are on the order of 1 meV. 




(10) 
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IV. RESULTS AND COMPARISON TO EXPERIMENTS 



A. The Quantum Wire 

1. Physical system 

The quantum wire considered here is a long, row-shaped, layered structure of fixed total height 
h. The geometry of the structure and the thickness and composition of each layer, based on the 
experimental work of Lukey et al. ^, are shown in Figure 7. The middle layers of the device 
are considered to be the active region, and include the quantum well layer (Sio.78Geo.22) and 
the two undoped barrier layers (Si). Surrounding the active region are strained (Sio.78Geo.22) 
emitter and collector regions; the thick outermost layers of the device are Si. A range of widths 
w arc considered in order to model size-dependent strain effects and to compare results with 
experimental data. 

The device operates by the sequential tunneling of charge carriers from the doped region above 
the barrier layer, into the quantum well layer, and then through the lower barrier layer. Resonant 
tunneling spectroscopy is done experimentally by applying a bias across the device and measuring 

the current of tunneling carriers that is induced. A resonant tunneling spectrum can be compiled 
by measuring the induced currents associated with a range of applied biases. The experimental 
result is a resonant tunneling I(V) curve. 

2. Strain field 

The strain and displacement fields for this geometry are two dimensional since the constraining 
effect of the material in the direction along the long axis of the wire imposes a state of plane 
strain. The mesh is more refined in the active region of the structure and near the traction free 
lateral surfaces, where the deformation is expected to be more nonuniform. The €22 component 
of strain (extensional strain in the X2 direction) is shown in Figure 3 in and near the active region 
of the device, plotted on a deformed mesh that shows the actual displacements, magnified for 
clarity. In Figure 4 the en component of strain along the center line of the quantum well layer 
is plotted for structures with three different widths. The important features are that the strain 
is a tensor valued function and that the relaxing effect of the free surface and the multilayered 
composition of the structure lead to nonuniform strain components. 

3. Results of the quantum mechanical calculation 

Solution of the quantum mechanical problem gives the energies and wavefunctions of states con- 
fined to the quantum well layer in the wire. From this spectrum of states, a density of confined 
states is calculated. The density of confined states for a narrow wire (w=250nm) is shown in 
Figure 8, representative eigenstates are shown in Figures 9 and 10, and the density of confined 
states for wide and narrow wires are compared in Figure 11. 

The results of the calculations are consistent with some experimental observations of Lukey 
et al. . First, the strain separates the resonant current peaks associated with the heavy hole and 
light hole subbands, as shown in Figure 8. The calculation gives a peak separation of about 20 
meV, which would correspond to a bias shift of roughly 50 mV. The experimentally measured 
separation is approximately 90 mV. However, it is important to note that the calculation shown 
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is for a wire aligned along a (100) crystalline axis. The experiments of Lukey et al. measured 
wires aligned along a (110) axis, which would exhibit more sensitivity to strain in the electronic 
properties due to the form of the deformation potential tensor D'^^ , and thus presumably a wider 
strain induced peak separation. The second characteristic evident in Figure 8 is the presence of 
additional fine features in the density of states. Fine structure is also observed experimentally in 
the w=250nm device. An examination of the states present over a range of energies shows that 
the fine features in the density of states are due to groups of similar states separated in energy 
by the influence of nonuniform strain, as shown in Figures 9 and 10. 

The size dependence of the strain effect is demonstrated in the densities of confined states 
for wires with widths of 250nm and 900nm in Figure 11. The main feature is the increase in the 
energy separation of the heavy hole and light hole peaks as the wire width increases. This is due 
to the reduced effect of free surface strain relaxation in larger devices, where the average strain 
values approach the bulk film mismatch strains. In the smallest devices, the strain is relaxed over 
a significant portion of the volume, so the average strain is smaller and the energy separation 
between the heavy hole and light hole peaks is smaller. 

B. The Quantum Dot 

1. Physical system 

The quantum dot considered here is cylindrical in shape and of fixed height. The features of this 
structure are based on the experimental work of Zaslavsky et al. ^ and Akyiiz et al. ^ and are 
shown in Figure 12 for a representative calculation. The three middle layers, which consist of 
the quantum well (Sio.75Geo.25) and the barriers (Si), make up the active region of the device. 
Surrounding the active region are the emitter and collector regions (Sio.75Geo.25) which have 
relatively low strain due to the outermost layers which are graded in composition. A range of 
cylinder diameters d is considered. 

The quantum dot device operates on the same resonant tunneling principle as the quantum 
wire. Carriers tunnel sequentially from the emitter region, through the upper barrier into states in 
the quantum well, and then through the lower barrier into the collector region. The experimental 
I(V) curve is a measure of the resonant tunneling spectrum. 

2. Strain field 

The stress, strain, and displacement fields are axisymmetric for this geometry. The mesh used 
to calculate the strain extends from the center axis of the structure outward, with increasing 
refinement near the outer, traction free surface of the device, where the fields are expected to be 
more nonuniform. Figure 13 shows the 622 component of strain (extensional strain in the axial 
direction) on a deformed mesh, and Figure 14 shows en, the extensional strain along a radial 
line in the midplane of the quantum well. As in the case of the quantum wire, the strain is very 
nonuniform, and the extent of the nonuniformity increases in smaller structures. This is due to 
the more significant effect of free surface proximity. 
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3. Results of the quantum mechanical calculation 

For the three-dimensional quantum dot, a reduced quantum mechanical basis is adopted in order 
to limit the total number of degrees of freedom in the calculation. To model only the lowest energy 
heavy hole states, it is sufHcient to consider a one-dimensional quantum mechanical basis where 
only the ||,+| > band is examined, but it is still possible to consider an anisotropic effective 
mass. 

Two significant results compare favorably to the experiments of Zaslavksy et al. ^ and Akyiiz 
et al. ^. First, as seen in Figure 15, there is a size dependence of the strain effect on the resonant 
energies that is similar to the effect in the quantum wires. In larger diameter devices, the average 
strain is higher because free surface relaxation is less significant, so the peak is shifted further 
from an idealized case without strain. 

Second, in the smallest devices, there are additional fine features in the density of states that 
are shown in Figure 16 to be the result of the influence of nonuniform strain. In agreement with 
the suggestion of Akyiiz et al. ^, it is evident that the relaxed strain near the lateral surface 
leads to low energy ring-like heavy hole states. The strain induced energy shifting of groups of 
states produces features in the density of confined states that are consistent with the I(V) curves 
for devices of the same size. Figure 17 shows a calculated density of confined states plot and a 
measured I(V) curve for a d=250nm quantum dot with lOmeV on the energy axis equal to 25 mV 
on the bias axis. Many of the features of the I(V) curve arc predicted qualitatively in the density 
of confined states curve, including fine structure below and above the main heavy hole resonance 
energy. However, the density of confined states is not equivalent to the resonant tunneling current; 
the calculation does not consider some important physical effects, most notably the roughly linear 
background current in the I(V) relationship. 

V. CONCLUSIONS 

A finite element technique is presented here which allows for the calculation of strain effects on 
the electronic and transport properties of strained quantum wires and dots. The approach is 
similar to some recent work as it is based on a simplified quantum mechanical model^"^, but 
the fiexibility of meshing and the low computational cost of the finite element method offer easy 
access to results which can be compared to experimental data. 

The technique is used to examine mismatch strain effects in quantum wires and quantum dots 
that operate on a simple single carrier sequential tunneling effect. Strain effects are shown to 
explain several reported trends in experimental data. In particular, two features of the mismatch 
strain in the devices have strong effects on the calculated electronic and transport properties. 
First, the average effect of the strain is to separate the resonant energy peaks associated with the 
individual valence subbands in the material. In larger devices, the strain is less relaxed by the 
free surfaces, so the HH-LH subband energy separation is larger. Second, strain nonuniformity in 
the devices is responsible for fine structure in the resonant tunneling current peaks. This effect 
is the source of low energy ring-like states that are found in the small cylindrical quantum dots, 
and the edge states found in quantum wires over a range of sizes. 

The main weaknesses in the method are in the simplified quantum mechanical model. The 
real-space calculation is based on a A;-space material model that is accurate near k = Q. The 



10 



resonant tunneling model assumes ballistic transport of a single charge carrier, and contact is made 
with experiments in only an approximate way. A linear elastic constitutive model provides a good 
approximation for the strain, although the approach has limitations for such small, highly strained 
structures'^. Finally, the strain effect is treated as a linear perturbation to a perfect crystal 
Hamiltonian, so the fully coupled nature of strain and electronic properties through the chemical 
bonding is not considered. However, the technique shown here is a promising, computationally 
inexpensive way to determine strain effects on electronic properties in semiconductors. The means 
to overcome the noted shortcomings are under development. 
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Appendix A 



The strain induced potential V^"^^(r) is given by 
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where each component D'^^{r) of the matrix for fixed forms a scalar product with the strain 
tensor eij{r) through summation over i and j. And similarly, the k • p Hamiltonian given by 
Luttinger and Kohn takes the form 
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where the each of the matrix components L'^^r) for fixed a/3 form a scalar product with the oper- 
ator V^^-.The components D'^^if) and Ll^^if) have very similar form. The deformation potential 



components Dff^f) are 
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D]:f{f) = Df^(f) = D^(f) = Dtf(f) 
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The Hamiltonian components L°'-^{r) can be obtained by making the substitutions 2m^'^i^^ 

ij 



72<->6, 73<-^<i into the expressions for the components D'^f{f), where 71, 72, 73 are the 



mo 
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Luttingcr-Kohn parameters. Values for the deformation potential constants and the Luttinger- 
Kohn parameters for Si and Ge are given below. Values for alloys of Si and Ge are interpolated 
from values for the bulk materials. 





a{eV) 


b{eV) 


d{eV) 


71 


72 


73 


Si 


2.1 


-1.5 


-3.4 


4.29 


0.34 


1.45 


Ge 


2.0 


-2.2 


-4.4 


13.4 


4.24 


5.59 



Appendix B 



To obtain the finite element form of the Schrodinger equation, the physical region is divided into 
elements, which are taken here to be 4-nodcd quadrilaterals for the two dimensional quantum wire 
and 8-noded bricks for the three dimensional quantum dot. The wavefunction wavefunction 
gradient V*°, and potential V^^ are expressed in terms of discretized values at the nodes, and 
values within the elements are determined by linear interpolation using linear shape functions 
N{r). A group of quadrilateral elements and the linear shape function for an associated node is 
shown in Figure 18. 

The form of the Schrodinger equation to be solved is 



{Bl) 



The weak form of the equation is obtained by multiplying by ^(r) and integrating over the volume 
of the body. The first term is integrated by parts, and the functional corresponding to the weak 
form is given by 



{B2) 



Jr Jr Jr 

The spatially varying fields are then discretized using the shape functions N{r) to get 
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m 

The total variation of the functional n(*^) is then minimized with respect to the nodal values 
of the wavefunction ^' □ so that 



A B 



thus 
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Replacing integrals over the region R with integration over individual element volumes (ri*^) and 
a summation over all elements, the final finite element form of the equation becomes 
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The contribution of a single element to the left hand side of the equation is given by the clement 
stiffness matrix. The integration over the element is done numerically at a set of quadrature 
points. The construction of an element stiffness matrix for the case of two spatial dimensions and 
a four subband quantum mechanical basis is as follows: 

= j^^ (vNaL^^VNb + NaNb ^ vfNc^ dn 

int.pts. nodes 



16x8 8x8 8x16 16x4' 7"';; '4x16 

4x4 



where the tensors D'^^ and L'^^ given in Appendix 1 reduce to 2 x 2 matrices, and the shape 
function matrix is given by 
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and the shape function derivative matrix VA^^ follows in the same form. 

The contribution of a single element to the right hand side of equation (B6) is referred to 
as the element mass matrix and is constructed in a similar way as k^^. The final finite element 
matrix form of the Schrodinger equation, given by 



is constructed by assembling the clement stiffness matrices and element mass matrices into global 
element and mass matrices, Kij and Mjj, for the entire device. 
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List of Captions 



1. Schematics of the quantum wire and quantum dot. The row-shaped quantum wire (Lukey 
et al. ) has width w and extends a distance much larger than w in the (010) direction. The 
cyhndrical quantum dot (Akyiiz et al. ) has diameter d. 

2. Finite element mesh used to calculate strain in the quantum wire. Since the strain is 
symmetric laterally and vertically about the center axes of the device, it is possible to 
model one quarter of the cross section of the structure only. A portion of the highly refined 
mesh near the edge of the active region is shown. 

3. Direct vertical component of strain in the upper half of the quantum wire, from the center 
to the edge. Displacements of the free outer surface of the device can be seen at the right 
edge of the plot. The layers, from the bottom, represent the well and barrier layers (highly 
nonuniformly strained) and the emitter and topmost Si layers. 

4. Direct lateral component of strain in the center of the well layer of the quantum wire. The 
strain is uniform near the center of the wire (left) but highly nonuniform near the edge 
(right). In larger wires the strain is less relaxed and less nonuniform. 

5. The total potential V{r) for the heavy hole band in the strained quantum dot. The potential 
in the barrier layers is high relative to the potential in the quantum well layer; the potential 
is axisymmctric, and radially nonuniform throughout the device. 

6. Representative probability densities in a d=50nm quantum dot. On the left is a low energy 
state with six-fold angular quantization localized in a ring-like region near the outer edge of 
the active region. On the right is a higher energy state, localized in the center of the active 
region, with two fold angular quantization and two fold quantization in the z direction. 

7. Schematic of the quantum wire geometry and composition in the active region. The 59A 
Si barriers surround the 33A Sio.78Geo.22 quantum well layer. An applied bias V induces a 
tunneling current in the z direction. 

8. Density of confined states in the strained quantum wire of width 250nm. The two large 
peaks in the dashed curve show the calculated heavy hole and light hole resonances without 
considering strain effects. Strain causes the resonances to separate in energy and induces 
fine structure in the density of confined states. States at points A (edge state) and B (light 
hole state) are shown in Figures 9 and 10. 

9. Probability density for a lower energy confined state. The four cross sections of the active 
region show the probability density associated with each of the four valence subbands for a 
given cigenstate. The wavefunction in this eigenstate is of mixed type, and localized near 
the edge of the structure due to the free surface strain relaxation. 

10. Probability density for a higher energy confined state. The predominantly > type 

wavefunction is localized in the center of the quantum well, but shows some effects of the 
nonuniform potential near the lateral surfaces and strong mixing with the ||, — | > subband. 
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11. Densities of confined states for two quantum wires of different widths. The higher average 
strain in the wider device results in a larger energy separation between heavy hole and light 
hole resonance peaks. The heavy hole peak is shown to be shifted by the strain more than 
the light hole peak. 

12. Geometry and composition of the quantum dot device (Akyiiz et al. ). The 45A Si barriers 
surround the 35A Sio.75Geo.25 well. The total height is 80nm, and the diameter d varies. 
An applied bias V induces a tunneling current in the z direction. 

13. Vertical component of strain in the upper half of the quantum dot, from the center to the 
edge. The plot is deformed to show the displacement in the structure, and the strain is 
highly nonuniform, particularly in the active region of the device. The two layers at the 
bottom are the quantum well and top barrier layers. 

14. Direct lateral component of strain in the well layer of the quantum dot, from the center to 
the edge. The strain is more relaxed and more nonuniform in the smaller structures. The 
strain profile is very similar to the profile in the quantum wire. 

15. Densities of heavy hole confined states for a range of quantum dot diameters. The variation 
in average strain levels results in the shifting of the peaks for dots of different sizes. The 
energy shift corresponds to the bias shift measured by Zaslavsky et al. . 

16. The effect of strain induced lateral confinement on the density of confined states. The 
features of the density of confined states are due to groups of eigenstates with similar lateral 
quantization. 

17. Fine structure in the density of confined states in the d=250nm quantum dot (Akyiiz et 
al. ). The calculated density of confined states is consistent with the experimental current- 
voltage curve. The density of confined states calculation does not account for the roughly 
linear background current that is observed experimentally. 

18. Elements, nodes, and a representative shape function in the two dimensional finite element 
formulation. 
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